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ABSTRACT 

We report the discovery of broad Wolf-Rayet emission hnes in the Multiple Mirror 
Telescope (MMT) spectrum of the NW component of I Zw 18, the lowest-metalhcity blue 
compact dwarf (BCD) galaxy known. Two broad WR bumps at the wavelengths A4650 
and A5800 are detected indicating the presence of WN and WC stars. The total numbers 
of WN and WC stars inferred from the luminosities of the broad He II A4686 and C IV 
A5808 lines are equal to 17±4 and 5±2, respectively. The WR-to-0 stars number ratio 
is equal to ~0.02, in satisfactory agreement with the value predicted by massive stellar 
evolution models with enhanced mass loss rates. The WC stars in the NW component 
of I Zw 18 can be responsible for the presence of the nebular He II A4686 emission line, 
however the observed intensity of this line is several times larger than model predictions 
and other sources of ionizing radiation at wavelengths shorter than 228A are necessary. 

Spectroscopic observations presented herein were obtained with the Multiple Mirror Telescope, a facility operated jointly by the 
Smithsonian Institution and the University of Arizona. 

Visiting astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatories, operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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1. INTRODUCTION 

The presence of large numbers of Wolf-Rayet (WR) 
stars in star-forming galaxies is well established (Allen, 
Wright & Goss 1976; D'Odorico & Rosa 1981; Oster- 
brock & Cohen 1982; D'Odorico. Rosa & Wampler 
1983; Hutsemekers & Surdej 1984; Kunth & Jou- 
bert 1985; Kunth & Schild 1986; Sargent & Filip- 
penko 1991; Conti 1991; Vacca & Conti 1992). These 
galaxies are often called WR-galaxies and they are 
of quite heterogeneous types. We focus here on 
the problem of detection of WR stars in the low- 
mass and low-metallicity blue compact dwarf (BCD) 
galaxies. Systematic spectroscopic studies of BCDs 
have shown that in the spectra of ~l/3 of BCDs, 
broad WR bumps characteristic of late WN stars are 
present, mainly at A4650 (Izotov, Thuan & Lipovet- 
sky 1994, 1997; Izotov & Thuan 1997a). The intensity 
of these bumps decreases with decreasing metallicity, 
in agreement with predictions of massive star evolu- 
tion models and models of evolutionary population 
synthesis for star-forming regions (Maeder & Meynet 
1994; Cervino & Mas-Hesse 1994; Leitherer & Heck- 
man 1995; Meynet 1995; Schaerer 1996). The lowest 
metallicity BCD in which WR stars have been de- 
tected has ~Zq/10, although BCDs can be as metal- 
deficient as Z0/5O. Are WR stars present in these 
extremely metal-deficient BCDs? In principle, mas- 
sive stellar evolution theory (e.g. Maeder fc Meynet 
1994) does predict the evolution of the most massive 
low-metallicity stars through the WR stage. How- 
ever, since the efficiency of mass loss by stellar wind 
decreases with decreasing metallicity, the total num- 
ber of WR stars and the total duration of the WR 
phase in a star formation episode are significantly re- 
duced at low metallicities. This trend led Schaerer et 
al. (1997) to conclude that in some metal-deficient 
BCDs, weak WR spectral features are not detected 
simply because of inadequate signal-to- noise ratio. 

Several recent observations of low-metallicity BCDs 
have suggested that massive stars with mass loss are 
indeed present in galaxies with heavy element abun- 
dances less than Zq/20. Imagery of I Zw 18 with 
HST by Hunter & Thronson (1995) has resolved its 
NW and SE components into stars, with the brightest 
star having ^ ^ 22 mag. Those authors attempted 
to find WR stars using a narrow-band image in the 
He II A4686 line and detected only two marginal 
WR candidates. They expected a large population 
of WR stars, given the presence of a large number of 



red supergiants, so concluded that I Zw 18 is not a 
BCD with many Wolf-Rayet stars. Recently, Izotov 
& Thuan (1997b), from 4m Kitt Peak Mayall tele- 
scope spectrophotometry of I Zw 18, noted that WC 
stars are possibly present in the NW component. Fi- 
nally, Thuan & Izotov (1997) have found evidence 
for stars with mass loss through the presence of P 
Cygni profiles in the UV HST spectra of two other 
very metal-deficient galaxies, SBS 0335-052 (^0/40) 
and Tol 1214-277 (^0/23). In this paper we continue 
our search for stellar populations with mass loss in 
very metal-deficient BCDs. We present high signal- 
to-noise ratio optical spectrophotometry of I Zw 18. 
We show that WR stars of different types are clearly 
present in this galaxy. 

I Zw 18 is a BCD undergoing an intense burst of 
star formation. It was first recognized to have an 
exceptionally low metal abundance by Searle & Sar- 
gent (1972). Later studies by Lequeux et al. (1979), 
French (1980), Kinman & Davidson (1981), Pagel 
et al. (1992), Skillman & Kennicutt (1993), Martin 
(1996), Izotov, Thuan & Lipovetsky (1997) and Izotov 
& Thuan (1997b) have confirmed the oxygen abun- 
dance to be only ~l/50 of the solar value. Zwicky 
(1966) described I Zw 18 as a double system of com- 
pact galaxies, which are in fact two bright centers 
of star formation separated by an angular distance 
of 5'.'8. These two star-forming regions will be re- 
ferred to subsequently as the brighter NW and fainter 
SE components. Later studies (Davidson, Kinman & 
Friedman 1989; Dufour & Hester 1990) have revealed 
a more complex structure with several additional dif- 
fuse features. However the NW and SE components 
dominate in brightness and this paper will focus on 
them. We describe the observations and data reduc- 
tion in §2. In §3 we discuss the properties of WR stars 
in the BCD. We summarize our findings in §4. 

2. OBSERVATIONS AND DATA REDUC- 
TION 

Spectrophotometric observations of I Zw 18 were 
obtained with the Multiple Mirror Telescope (MMT) 
on the nights of 1997 April 29 and 30. Observations 
were made with the blue channel of the MMT spec- 
trograph using a highly optimized Loral 3073x1024 
CCD detector. A l'.'5xl80" slit was used along with 
a 300 groove mm^^ grating in first order and an L- 
38 second-order blocking filter. This yields a spatial 
scale along the slit of 0'.'3 pixel^^, a scale perpendicu- 



3 



lar to the slit of 1.9 A pixel" ^, a spectral range of 3600 
- 7500A, and a spectral resolution of 7A (FWHM). 
For these observations, CCD rows were binned by a 
factor of 2, yielding a final sampling of C'G pixel"^. 
The observations cover the full spectral range in a 
single frame that contains all the lines of interest and 
have sufficient spectral resolution to distinguish be- 
tween narrow nebular and broad WR emission lines. 
The total exposure time was 180 minutes and was bro- 
ken up into six sub-exposures, 30 minutes each. All 
exposures were taken at small air masses (1.1 - 1.2), 
so no correction was made for atmospheric dispersion. 
The seeing was 0'.'7 FWHM. The slit was oriented in 
the direction with position angle P.A. = -41° to per- 
mit observations of both NW and SE components. 
The spectrophotometric standard stars EG 247 and 
HZ 44 were observed for flux calibration. Spectra of 
He-Ne-Ar comparison lamps were obtained after each 
subexposure to provide flux calibration. 

The two-dimensional spectra were bias subtracted 
and flat-field corrected using the IRAF ^. For the 
NW component, the extracted one-dimensional spec- 
tra cover the brightest part of the galaxy with a spa- 
tial size of ~5". Similar procedures were used for the 
SE component, resulting in one-dimensional spectra 
covering a region 5" wide at a distance of 5'.' 8 from 
the NW component. The extracted spectra from each 
frame were then coadded and calibrated to absolute 
fluxes. 

The observed line intensities have been corrected 
for interstellar extinction using the reddening law 
by Whitford (1958). Hydrogen lines have been also 
corrected for underlying stellar absorption, with the 
equivalent width for hydrogen absorption lines de- 
rived self-consistently together with the extinction co- 
efhcient from the observed intensities of all hydrogen 
lines. We show the observed F{X)/ F(ll(3) and ex- 
tinction and absorption-corrected /(A)//(H/3) line in- 
tensities for the NW and SE components in Table 1, 
along with the extinction coefficient (C((H/3)) and the 
equivalent width of the hydrogen absorption lines, the 
observed flux and equivalent width of the H/3 emis- 
sion line. The uncertainties for the tabulated relative 
line intensities are ~0.5% for the strongest lines and 
'^lO - 15% for the weakest lines. 



•^IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 



3. WOLF-RAYET STARS IN I ZW 18 

In Figure 1 we show the spectra of the NW and the 
SE components of I Zw 18. The WR broad lines and 
strong narrow He II A4686 emission line are clearly 
present in the spectrum of the NW component. Thus, 
I Zw 18 is the lowest metallicity galaxy where WR 
stars are detected. In contrast, in the spectrum of 
the SE component WR emission lines are not con- 
vincingly detected. The weak emission seen at the 
wavelength ~ 4700 A is nebular emission from the 
[Fe HI] A4658, He II A4686 and [Ar IV] AA4711, 4740 
lines. The appearance of a broad underlying emission 
feature here most probably arises from the blending 
of the wings of the individual nebular line profiles, 
but a very small contribution from WR stars can- 
not be excluded. The detection of the C IV A5808 
broad line in the NW component suggests the pres- 
ence of early WC stars, while the other broad lines 
are, most likely, evidence of early and late WN stars. 
The formation of C IV A5808 by WN stars can be 
ruled out because of its large FWHM ('^80A) which 
corresponds to early type WC stars (Smith, Shara 
& Moffat 1990). The approximate number of O and 
WR stars in the NW component can be estimated as 
follows. The observed fiuxes of the C IV A5808 and 
broad He II A4686 emission lines are 8.3x10^^^ and 
1.49x10-1^ erg s'^cm'^, and that of H/3 A4861 inte- 
grated along the slit is 4.93x10"^"* erg s~^cm~^. Cor- 
recting for interstellar extinction C(H/3) = 0.13 dex 
for the NW component and adopting a distance D = 
10.8 Mpc (for the redshift 2^=0.00274 of the NW com- 
ponent with Ho = 75 km s^^ Mpc~^), we derive the 
following luminosities: L{C IV A5808) = 1.27x10^'^ 
erg s-\ L(He II A4686) = 2.84x10^^ erg s^^ and 
L(H/3) = 9.26x10^* erg 3"^ Assuming that only 1/2 
of the light in H/3 is contained within the slit width 
of I'.'S, we derive finally L(H/3) = 1.85x10^9 erg s'^ 
This value agrees well with that derived by Hunter & 
Thronson (1995) from the HST Ha image. Adopting 
a value of Qo = 49.05 (Vacca & Conti 1992) for the 
logarithm of the number of Lyman continuum pho- 
tons emitted per second by an 07V star and assuming 
Case B recombination and an instantaneous burst of 
star formation, we derive A^(07V) — 381. To derive 
the total number of O stars, we need to take into ac- 
count the age of the stellar population and the IMF 
slope. A measure of the age of the burst of star for- 
mation is the equivalent width of the H/3 emission 
line. The low equivalent width of 56A in the NW 
component of I Zw 18 corresponds to an age of 4 - 5 
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Myr, assuming a Salpctcr IMF (Lcntlic'rc;r & Hcckman 

1995) . This value is in good agreement with the re- 
sults of direct HST photometry of the brightest stars 
by Hunter & Thronson (1995). They fomid stars as 
old as 5 Myr in the NW component. Then, the to- 
tal number of O stars is ~ 3 times larger (Schaerer 

1996) and is equal to ^1100. This approximate value 
is in good agreement with the number of O stars of 
1300 derived by Hunter & Thronson (1995) from the 
luminosity function. 

The number of WR stars is estimated from the 
luminosity of the WR emission features. While the 
main contribution to the A5808 emission is from 
WC4 stars, the emission in the A4650 feature is pro- 
duced by WR stars of different types. However, the 
observed flux ratio, f(A4650)/f(A5808)-3, is signif- 
icantly greater than that expected for WC4 stars 
(Smith, Shara, & Moffat 1990). Therefore, the domi- 
nant contributors to the broad He II A4686 emission 
are WNL stars. The moderate spectral resolution and 
line blending, however, prevent drawing more decisive 
conclusions about the origin of the A4650 bump. The 
derived number of WNL stars should therefore be con- 
sidered as only indicative. Adopting the luminosity 
of a single WC4 star in the C IV A5808 line as equal 
to 2.5x10^6 erg s-\ and that of a single WNL star in 
He II A4686 as equal to 1.7x 10^6 erg s'^ (Conti 1991; 
Vacca & Conti 1992), wc derive the following numbers 
of stars: 7V(WC4) = 5±2, iV(WNL) = 17±4. This 
gives 7V(WR)/7V(0) = 0.02 and 7V(WC)/7V(WN) = 
0.3. 

Meynet (1995) presented new evolutionary popula- 
tion synthesis models based on the most recent grids 
of stellar models computed at the Geneva Observa- 
tory. He studied the effects of changes in the rates of 
mass loss by stellar winds on the massive star popu- 
lations born in a starburst. According to Maeder & 
Meynet (1994), the high mass loss rate stellar models 
are to be preferred over the standard ones on the basis 
of comparisons with the observed luminosities, chem- 
ical compositions and nmiibcr statistics of WR stars 
in zones of constant star formation rate. In starburst 
galaxies, the presence of WC stars at very low metal- 
licity is predicted only by models of massive star evo- 
lution with enhanced mass loss (Meynet 1995), while 
the models with standard mass loss rates derived by 
de Jager, Nieuwenhuijzen & van der Hucht (1988) and 
scaled with metallicity as fail to produce WC 
stars at the metallicity of I Zw 18. Therefore, the 
detection of WC stars in I Zw 18 gives strong sup- 



port to the idea of enhanced mass loss in massive 
low metallicity stars. Furthermore, we find satisfac- 
tory agreement between the observed and theoreti- 
cal WR/0 and WC/WN ratios at the metaUicity of 

I Zw 18. Meynet (1995) has calculated evolution- 
ary population synthesis models only for metallici- 
ties as low as .^0/20. At this metallicity, assuming 
an instantaneous burst of star formation and an IMF 
dN/dMocM~'^ for the massive stars, his models pre- 
dict WR/OwG.035 and WC/WN«0.5 for peak values. 
Scaling these values as Z^-^ to the I Zw 18 heavy el- 
ement abundance Zq/50 gives WR/0 and WC/WN 
ratios close to those derived from the observations. 

The detection of WC stars in I Zw 18 can help re- 
solve the long-standing problem of the origin of the 
strong nebidar He II A4686 line in the NW compo- 
nent, which is several orders of magnitude greater 
than predicted by photoionized H II region models. 
Schaerer (1996) has shown that WC stars can signif- 
icantly increase the ionizing flux shortward of 228A , 
thus leading to the formation of a He++ zone in the 
H II region and increasing the recombination He II 
A4686 emission line luminosity by several orders of 
magnitude. However, this model predicts the maxi- 
mum value of the nebular He II A4686 emission line 
intensity to be only ^1% 2% that of H/3 when the 
age of the star forming region is not greater than 3 
Myr. The He II line intensity is expected to decrease 
from this maximum value for a somewhat older stel- 
lar population like the one in I Zw 18. This is in 
contrast to the observed intensity of the nebular He 

II A4686 emission line in the NW component of I Zw 
18 which is 4% that of H/3, larger than the predicted 
peak value. Furthermore, this model fails to explain 
the presence of strong He II A4686 in the galaxies 
where WC stars have not been detected (Izotov et 
al. 1997). We conclude therefore that, in addition to 
WC stars, some other source of ionizing radiation at 
wavelengths shorter than 228A must be invoked. 

4. CONCLUSIONS 

We have obtained high signal-to-noise ratio spec- 
trophotometric observations of I Zw 18, the most 
metal-deficient BCD known, in an attempt to detect 
the broad low-intensity emission lines of WR stars. 
We have obtained the following results: 

1. Broad emission lines at wavelengths ~ A4650 

and A5800 have been detected in the spectrum of the 
NW component of I Zw 18 implying the presence of 
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WN and WC stars, while in the younger SE com- 
ponent these Unes are absent. Thus, I Zw 18 is the 
lowest-metalUcity Wolf-Rayet galaxy known to-date. 

2. The total numbers of WNL and WC4 stars in I 
Zw 18 are 17±4 and 5±2 respectively, and the total 
number of O stars is ~ 1100. The existence of WC 
stars at the very low metallicity of Zq/50 in I Zw 
18 confirms predictions of massive stellar evolution 
models with enhanced mass loss rates (Meynet 1995), 
while the models with standard mass loss rates fail 
to explain the presence of WC stars in I Zw 18. The 
observed WR/0 and WC/WN ratios of 0.02 and 0.3 
respectively, are in satisfactory agreement with evolu- 
tionary population synthesis models based on stellar 
evolution models with enhanced mass loss and extrap- 
olatcxl to the; metallicity of I Zw 18. Although the WC 
stars detected in I Zw 18 could be responsible for the 
presence of the strong nebular He II A4686 emission 
line, the observed value of the nebular He II A4686, 
4% that of H/3, is several times larger than theoretical 
predictions. We conclude that additional sources of 
ionizing radiation at wavelengths shorter than 228A 
need to be present. 
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Fig. 1. The spectra of the NW and the SE compo- 
nents in the BCD I Zw 18. The broad nitrogen and 
carbon lines in the spectrum of the NW component 
are marked indicating the presence of WN and WC 
stars. The permitted Mg 1 A4571 and narrow nebu- 
lar He 11 A4686 emission lines are also marked. The 
broad WR emission Uncs arc absent in the spectrum 
of the SE component where the nebular [Fe 111] A4658, 
He II A4686 and [Ar IV] AA4711, 4740 emission Unes 
arc marked. The spectra have been smoothed using 
a 3-point box-car. 



Table 1 
Emission Line Intensities 



I Zw 18 (NW) 



I Zw 18 (SE) 



Ion 


F(X)/F(B.d) 








3727 [0 II] 


0.228 


0.238 


0.510 


0.499 


3868 [Ne III] 


0.143 


0.148 


0.139 


0.136 


3889 He I + H8 


0.102 


0.195 


0.155 


0.200 


3968 [Nc III] + H7 


0.134 


0.222 


0.169 


0.212 


4101 RS 


0.196 


0.276 


0.232 


0.270 


4340 H7 


0.410 


0.472 


0.443 


0.468 


4363 [0 III] 


0.069 


0.068 


0.054 


0.052 


4471 He I 


0.022 


0.021 


0.034 


0.033 


4686 He II (neb) 


0.041 


0.040 


0.007 


0.006 


4686 He II (WR) 


0.063 


0.062 






4861 H/3 


1.000 


1.000 


1.000 


1.000 


4959 [0 HI] 


0.731 


0.690 


0.591 


0.573 


5007 [0 HI] 


2.197 


2.069 


1.776 


1.724 


5808 C IV (WR) 


0.035 


0.031 






5876 He I 


0.074 


0.066 


0.095 


0.092 


6563 Ha 


3.163 


2.743 


2.834 


2.745 


6678 He I 


0.030 


0.026 


0.028 


0.027 


6717 [S II] 


0.024 


0.021 


0.043 


0.042 


6731 [S II] 


0.019 


0.016 


0.031 


0.030 


7065 He I 


0.026 


0.022 


0.024 


0.024 


7135 [Ar HI] 


0.019 


0.016 


0.019 


0.018 



C(H/3) dex 
F(H/3)^ 
EW(Rf3) A 
Ewlabs) A 



0.130 
2.36 
56 
2.9 



0.010 
1.76 
129 
3.9 



'^in units of 10 ergs s ^cm ^ 
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Wavelength (A) 
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